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Cálculo simplificado de parámetros ortodrómicos para navegación: 
Derrota ortodrómica y rumbo inicial. 


O 2022 Carlos Grasa Lambea 


En este documento deduciremos algunas fórmulas básicas para el cálculo en la navegación ortodrómica, 
pensando en usarlas en algunas aplicaciones para móviles. Se intentará la máxima simplicidad cubriendo unos 
mínimos preestablecidos de precisión en los resultados. El documento pretende ser didáctico mediante 
explicaciones e ilustraciones, a la vez que riguroso en la deducción de las fórmulas. Se incluirán conceptos 
básicos necesarios para el fácil seguimiento de la lectura. 


Sección 1: Cálculo simplificado de la distancia ortodrómica entre dos puntos de la superficie terrestre. 


Varias son las formulaciones existentes para calcular la distancia mínima entre dos puntos de la superficie 
terrestre sin abandonar ésta y sin desviarse a derecha o a izquierda, es decir: siguiendo una trayectoria 
geodésica. Si el lector está interesado en ojear otros desarrollos, los encontrará fácilmente en los temas: 
triángulo esférico, fórmulas de Bessel, de Haversine,... 


Aquí realizaremos un desarrollo básico a partir de conceptos fundamentales, como si se enfrentara este 
problema por primera vez. 


Para calcular la derrota ortodrómica, o distancia ortodrómica, calcularemos el arco de circunferencia 
comprendido entre los dos puntos de la superficie terrestre. Este círculo es llamado máximo porque su centro 
coincide con el centro de la tierra y es el mayor círculo que la esfera terrestre puede contener. 


Aquí estamos haciendo una primera aproximación al considerar la tierra como una esfera. Esto simplificará 
mucho los cálculos, pero tendremos que cuantificar el error que estamos dispuestos a asumir en el resultado. La 
tierra es más parecida a un elipsoide de revolución con un cierto achatamiento de los polos. Para más 
información, mírese el tema del geoide terrestre. 


Podemos visualizar nuestro objetivo observando el 
esquema del círculo de la ilustración. 


Queremos conocer la distancia d de la trayectoria del 
arco de círculo máximo que separa los puntos A y B 
cuyas posiciones son conocidas. Este círculo máximo 
lo es porque su centro coincide con el centro O de la 
esfera terrestre. Este origen O dista de los puntos A y 
B precisamente el radio R terrestre. Se nos ocurre 
calcular la longitud del arco d como el producto del 
radio R por el ángulo w formado por los radios hacia 
los puntos A y B. 


Así pues, obtenemos una primera fórmula, 


fórmula 1: d=R:'0w 
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El radio R es el radio terrestre, conocido, y tomaremos su valor promedio asumiendo perder un poco de 
precisión. El radio ecuatorial es de 6378 km, mientras que el radio polar es de 6356 km. 

Nosotros vamos a usar el radio volumétrico medio de 6371 km, que es el radio que tendría una esfera de igual 
volumen que nuestro planeta. Tomaremos R como 6371000 metros para nuestros cálculos. 


Pero ahora surge un objetivo intermedio porque no conocemos el ángulo w que forman los radios. Si volvemos 
a mirar la figura 1, observamos que los puntos A , B forman con el centro O un triángulo. Si conociésemos sus 
tres lados podríamos calcular el ángulo fácilmente con la fórmula para hallar el lado a de un triángulo 
conociendo los otros dos lados b y c, y el ángulo A entre ellos: 


c*=a?+b?-2:b:c-cos(A) 


En nuestro caso, c es la cuerda C, mientras que b y c son ambos el radio terrestre R, y el ángulo A es el que en la 
figura hemos llamado w en la figura 1. Y sustituyendo así obtenemos la fórmula 


C'=R+R*-2-R-R-c0sw y simplificando C'=2-R”(1—cos 0») 


Ignoramos la longitud de la cuerda C y el valor del ángulo w, pero si llegamos a calcular C podríamos obtener 
w despejándola y expresando la 


, G 
fórmula 2: w=arccos(1— 
(1-3) 


Ahora sustituimos w en la fórmula 1 por su expresión de la 
fórmula 2, y así conseguimos la 


fórmula 3: d=R:-arccos(1— 7 ) 


Nuestro problema consiste ahora en conseguir hallar la longitud de la cuerda C que es más sencillo, pues 
consiste en calcular la distancia recta entre dos puntos A y B del espacio. 


Si conociésemos las coordenadas cartesianas x, y, z de ambos puntos será inmediato el cálculo de la distancia 
que los separa. 


Nombramos las coordenadas cartesianas del punto A como: xa, ya, Za Correspondientes a los ejes cartesianos, e 
igualmente denominamos Xp, Ys, Za como las coordenadas del punto B. Ahora vemos claro que la longitud de la 
cuerda C que va desde el punto A hasta el punto B será el resultado de la 


fórmula 4: CO=[(x Xp) +H(Ya—Y 5) +(2A—=2p) 


Sólo un detalle: los puntos de la superficie terrestre suelen estar expresados mediante las coordenadas 
geodésicas de latitud y longitud. La latitud es el ángulo que hay que recorrer desde el ecuador terrestre hasta el 
punto de la superficie considerado. Todos los puntos del ecuador tienen latitud 0*, el polo norte tiene latitud de 
90” (positivos) y el polo sur tiene latitud de -90* (negativos). 


Cualquier punto del hemisferio norte tiene una latitud comprendida entre 0* y 90* (+90), mientras que los 
puntos del hemisferio sur tienen una latitud negativa entre 0* y -90* sexagesimales. 
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Hay una forma de expresar latitudes hacia el sur como, por ejemplo 45 %S para indicar 45 grados hacia el sur, 
pero vamos a expresar latitudes sur como si fueran latitudes norte negativas: -45 *N equivalentes. 

Por otra parte, la longitud es la otra coordenada geodésica que expresa el ángulo que hay que recorrer sobre el 
ecuador desde un meridiano de referencia (meridiano de Greenwich) hacia el este hasta llegar al punto del 
ecuador más próximo al punto de la superficie que deseamos referir. 


Hay varias maneras de expresar la misma longitud. Por ejemplo, 250* (hacia el este) sería equivalente a 110 %0, 
110” hacia el oeste. Pero nosotros usaremos siempre grados hacia el este. Cuando el ángulo sea mayor de 180 “E 
usaremos cantidades negativas hacia el este equivalentes a positivas hacia el oeste. Es decir, 110 “O será 
equivalente a 250 “E como también equivaldrá a -110 “E (grados este negativos). Nosotros usaremos esta última 
forma por ser la más actual y dominante. 


Generalmente la latitud y la longitud se expresan en grados sexagesimales. Aunque también hay otras unidades 
como grados centesimales o radianes. 


Estos últimos (radianes) son empleados por muchas calculadoras en sus funciones trigonométricas, por lo que a 
veces tendremos que convertir entre grados sexagesimales y radianes. Una vuelta completa tiene 360* pero 
también 2-1t radianes, por lo que 1? = 1 / 180 radianes. Y, a la inversa, 1 radián = 180 / rt grados. 


Tradicionalmente una fración de grado es expresada mediante minutos y segundos. De manera que 60 segundos 
hacen un minuto, y 60 minutos hacen un grado (aquí no una hora). El símbolo de grado es ” mientras que el de 
minuto es ' (comilla simple) y el de segundo es " (doble comilla). 


El uso de GPS y otros sistemas actualizados ha relegado el uso de los minutos y segundos para utilizar 
mayoritariamente los grados con cifras decimales. En este documento haremos uso del punto decimal (y no de 
la coma) para referirnos a números con parte fraccionaria. 


Convertir a esta notación es sencillo. Por ejemplo, supongamos que nos dan un ángulo (latitud o longitud) de la 
forma 24” 8' 35" y deseamos convertirlo a grados con decimales. Calcularemos 24 + 8 / 60 + 35 / 3600 para 
obtener el resultado 24.143055 grados sexagesimales, o también 24.143055* en forma más breve. 


Como detalle histórico aportamos que una milla naútica es la distancia recorrida sobre un meridiano cubriendo 
un ángulo de un minuto de latitud. Obsérvese que la distancia de un grado de longitud sobre un paralelo es 
indeterminada porque el radio de giro depende de la distancia al eje de la tierra y esta distancia depende a su 
vez de la latitud. 


Imagine que se halla a un metro del polo norte geográfico y hace un viaje de 180” hacia el este; la distancia 
recorrida será de tan solo de unos 3 metros y 14 centímetros, exactamente rt metros. 


Proseguimos con nuestra tarea de calcular el valor de la cuerda al cuadrado, es decir C? mediante la 
fórmula 4: C=lxa—xXa) + yr y a) +(Za—Zp) 


Antes tendremos que hallar las coordenadas cartesianas a partir de las respectivas latitudes y longitudes. Para 
ello vamos establecer qué símbolo representará cada uno de esos parámetros. Usaremos las letras minúsculas x, 
y, Z para expresar coordenadas cartesianas. Usaremos la letra 6 (delta) del alfabeto griego para significar latitud. 
Con la letra A (lambda) del alfabeto griego expresaremos la longitud. Y usaremos las mayúsculas 4 y ¿ como 
subíndices para expresar que un parámetro corresponde al punto A o al punto B. 
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La notación elegida queda así: 


Xa coordenada x del punto A 
Ya coordenada y del punto A 


Za coordenada z del punto A 
Xp coordenada x del punto B 
Y coordenada y del punto B 
Za coordenada z del punto B 


Óa latitud del punto A 
Ma longitud del punto A 
0 latitud del punto B 
Ag longitud del punto B 


Veamos ahora cómo convertir el sistema de 
coordenadas geodésicas (6, A) en un sistema de 
coordenadas cartesianas (x, y, z) para poder aplicar 
la fórmula 4 ya vista. Si nos fijamos en la figura 2, 
veremos que es muy fácil obtener la coordenada z 
del punto A usando de la trigonometría el seno de 
su latitud con la 


fórmula 5: z, =R-senó, 


Por otra parte el radio R proyecta sobre el plano 
ecuatorial una distancia horizontal H que luego 
usaremos para calcular otras coordenadas. Esta 
proyección H la calculamos con la 


fórmula 6: H=R:cosó, 


Ahora vamos a observar esta esfera desde una 
perspectiva cenital, situados encima del polo norte 
y mirando hacia el sur para observar el plano 
ecuatorial desde arriba. Esta situación queda 
representada en la figura 3. 


Recordemos que la circunferencia que vemos 
ahora es el ecuador terrestre. La proyección del 
radio al punto A sobre el ecuador es H y el ángulo 
corresponde a la longitud de A. Hemos escogido el 
eje X paralelo al meridiano O (de Greenwich). De 
nuevo mediante trigonometría calculamos las 
coordenadas Xa , Ya por medio de 


fórmula 7: x,=H-cosA, 


fórmula 8: y, =H:senA, 


eje Z 


eje Y 
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Añadimos una ilustración en 3D para 
redundar visualmente en lo pretendemos 
exponer. 


Intente comprender cómo obtenemos las 


coordenadas cartesianas Xa , Ya , Za del 
punto A en la figura 3b. 


Ahora en fórmula 7 y fórmula 8 
sustituimos H por la expresión de la 
fórmula 6 obteniendo 

fórmula 9: x,=R:c0sÓ,'CosA, 


fórmula 10: y, =R'cosó , send, 


A estas dos ecuaciones añadimos la ya 
conseguida 


fórmula 5: z, =R:senó, 
Y con estas tres fórmulas podremos pasar 
de latitud y longitud a coordenadas 


cartesianas. 


Ahora obramos análogamente con el punto B para obtener 


fórmula 11: Xg=R-:C0SÓ y' COSA 
fórmula 12: y¿=R:cos 0, send y 
fórmula 13: Z¿=R-senóz 


Estas seis fórmulas junto con las anteriores 


fórmula 4: Có=(xa— Xp) +(ya—Y nm) +(za—Z5) 


fórmula 3: d=R:-arccos(1— _ ) 


van a ser suficientes para calcular la distancia ortodrómica (o derrota ortodrómica) entre dos puntos A, B de la 
superficie terrestre, dada la rapidez con que los sistemas informáticos procesan el cálculo numérico. 


No obstante vamos a simplificar el cálculo intentando conseguir una única ecuación, más elegante. Si el lector 
no está interesado en esta siguiente sección, puede saltársela y quedarse únicamente con el resultado. 
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Pretendemos ahora conseguir una fórmula mucho más simple y directa para hallar la distancia ortodrómica 
entre dos puntos. 


Comenzamos reuniendo aquellas ecuaciones que serán la base de nuestro desarrollo: 


fórmula 4: Có=(xa=X8) +(Y a Y 5) +(24—2p) 
fórmula 9: x¿=R:c0sÓ ,'CosA, 

fórmula 10: y, =R:cosó , send, 

fórmula 5: z, =R-senó, 

fórmula 11: Xg=R-C0SÓ y' COSA y 

fórmula 12: y¿=R:cos 0, send y 

fórmula 13: z7¿=R:senóz 


Empezamos por sustituir las variables Xa , Ya, Za, Xp, Ys , Zp de la fórmula 4 por sus expresiones de las 
fórmulas 9, 10, 5, 11, 12 y 13: 


C*=(R-co0só,-cosA,—R-cosÓ ¿"cos Ag)" +(R-cos 8, send,—R-cosó¿"senAy)'+(R-send ,—R-sen Óg)' 


2 z .s . . . q. E 
Podemos sacar R” como factor común en la ecuación anterior y pasarlo dividiendo al otro lado de la misma, 
quedando: 


2 
ca=cos $,-c0sA,—c0s Óy"cos Ag) +(cos g,-senA,—cosÓ ¿sen Ag) +(R-senS ,—R-senóz)' 


A continuación vamos a desarrollar los cuadrados de los paréntesis del lado derecho de la ecuación: 


2 
C 2 2 2 2 
—G=c0s"Ó cos” A¿+COoS Óg:COS” Ag—2:C08Ó,*COSA,'COS Óg:COS A+ 
2 2 2 2 
+cos"Ó0 , sen” A,+cos" Ó sen” Ag—2:c0s Ó, send ,:COos Ó g:senAg+ 
2 2 
+sen* 0 ,+sen 9¿— 2:sen Ó ¿sen Óz 


Sacaremos factores comunes cos”g, , cos gg , COSÓ,:CosÓ, y quedará: 


2 
C 2 2 2 2 2 2 2 2 
paseos 8 "(cos A+sen”A,)+cos Ó(cos” Ag+sen"Ag)+sen*S +sen g¿— 


—2-c0s 9, cosÓ*(cosA cos Ay+senA,:senAg)—2-sen Ó ,: sen 5, 


Continuamos aplicando que la suma del seno al cuadrado más el coseno al cuadrado es la unidad. También 
aplicamos la fórmula del coseno de la diferencia: 


2 


C 2 2 2 2 
—G=co0s" 0 +cos 0g+sen' Ó ¿+sen $¿—2:c0s Ó,:cos 8: cos[1,—Az)—2:send ,senó y 
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Y seguimos simplificando: 


2 


=2-2:c08 5 ,c0sÓ ¿cos(A,—Ap)—2-sen 8, send 


Dividiendo ambos lados de la ecuación entre 2 y redistribuyendo los términos: 


2 
fórmula 14:  1- > Sen $, sen$y+c0só cos 9¿cos(1,—Ag) 


Si ahora recuperamos la ya deducida 


2 


) 


fórmula 3: d=R-arccos(1— 
2-R? 


Combinando ambas fórmulas concluimos con la expresión que nos permite obtener la distancia ortodrómica 
(derrota ortodrómica) entre dos puntos de la superficie terrestre, dadas sus coordenadas geográficas 
(geodésicas) de latitud y longitud: 


fórmula 15:  d=R:-arccos(sen S,:senóS¿+cosó,:cos Íy"cos(A4—Ag)) 


Mirando por encima la fórmula 15 apreciamos que no presenta problemas de estabilidad en cuanto a divisiones 
por cero o valores extremos. La expresión entre paréntesis (el argumento del arcocoseno) debe ser un número 
acotado desde -1 hasta +1 con ambos extremos incluidos, por lo que podríamos hacer una simple comprobación 
para limitar ese valor antes de calcular el arcocoseno. 


Esta fórmula es bien conocida y ampliamente usada para hallar la distancia ortodrómica. Nosotros la hemos 
deducido, invitando al lector a ser atrevido para realizar sus propios desarrollos aunque llegue a resultados que 
otros ya lograron, como en este caso. 


Por otra parte, la distancia calculada es directamente proporcional al valor escogido para el radio terrestre. 
Siendo el radio polar de 6356000 metros y el radio ecuatorial de 6378000 metros, hallamos una diferencia de 
22000 metros. Dividiendo esta diferencia entre el radio volumétrico medio que estamos usando de 6371000 
metros nos da 0.345314% como estimación del error relativo. Es decir, cada 100 km podemos cometer un error 
de 345 metros. Esto puede resultar aceptable para muchas aplicaciones que no necesiten posicionar con mucha 
exactitud objetos sobre el terreno, sino sólo informar sobre trayectos de naves en tiempo real. 


Concluimos aquí esta primera sección habiendo conseguido con la fórmula 15 un método sencillo para el 
cálculo de la distancia ortodrómica. 
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Sección 2: Cálculo simplificado del rumbo inicial ortodrómico entre dos puntos de la superficie terrestre. 


En esta nueva sección abordamos otro asunto básico de la navegación ortodrómica. Deseamos realizar una 
trayectoria, sin abandonar la superficie terrestre, lo más directa posible (ortodrómica) desde un punto inicial 
hasta otro punto final; ¿qué rumbo tenemos que establecer inicialmente al comienzo de nuestro viaje? 


El rumbo es la medida del ángulo en sentido horario tendido desde la dirección norte hasta la dirección en que 
avanzamos. Este parámetro es una coordenada horizontal que también llamamos azimut (o acimut). 


El azimut lo vamos a expresar en grados sexagesimales con decimales, siempre positivo, desde 0” hasta 360" 
donde 0* está incluido pero 360% no, porque es una redundancia de 0” a efectos de notación. 


Fijémonos ahora en la figura 4 para imaginar la situación. 
Estamos en el punto inicial y vamos a viajar por la 
trayectoria más corta hasta el punto final. Si conocemos el 
rumbo inicial Ri del viaje, primero apuntaríamos hacia el 
polo norte; luego rotaríamos en sentido horario tantos 
grados como el azimut dado, o rumbo inicial. Ahora 
avanzaríamos hacia nuestro destino. 


Observe que el azimut de llegada será generalmente 
distinto, así como todos los intermedios. Si realizamos la 
ruta ortodrómica, la más directa posible, veríamos variar 
la aguja de nuestra brújula contínuamente. 


En la práctica se establecen unos cuantos puntos 
intermedios consecutivos y solo se varía la dirección 
cuando alcanzamos la siguiente meta parcial. Entre dos 
puntos intermedios consecutivos no hay ajuste de 
dirección. 


Con la automatización de los sistemas de navegación, las naves pueden realizar decenas de ajustes de dirección 
cada minuto. Una aplicación sencilla puede determinar su nueva posición cada segundo mediante GPS y 
recalcular un nuevo rumbo inicial para el tramo de trayectoria que falta desde ese punto. 


Nosotros en esta sección vamos a deducir una fórmula que calcule el rumbo inicial de manera sencilla y que 
resulte práctica. De nuevo asumimos las imprecisiones de considerar la tierra como una esfera perfecta sin 
achatamientos. 


Hay otras fórmulas que se vienen usando desde hace mucho tiempo. Una de las más utilizadas recurre a la 
función arcotangente en el paso final del cálculo. Se trata de la conocida 


senlAz—Aa) 
tg(S,)-cos(S,)-senA,cos(8,—9,) 


) 


fórmula 16: — R¡=arctan( 


Esta formula presenta algunas limitaciones que ahora expongo. Se calcula la función arcotangente de una 
división cuyo denominador puede tender a cero, creando resultados que desborden la capacidad de los 
procesadores aritméticos. También vemos que en el primer término del denominador se calcula una función 
tangente cuyo resultado también puede tender a infinito. 
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Los problemas pueden aperecer cuando se da alguna de estas estas circunstancias: punto final cerca de alguno 
de los polos (la tangente tiende a +infinito); también hay ocasiones en que los puntos inicial y final se hallen en 
un mismo meridiano, produciendo una indeterminación 0 dividido entre 0. 


Si nos limitamos al cálculo de trayectorias no demasiado largas, menores de 10000 km, podemos deducir una 


fórmula que evite estos problemas. 


Vamos a considerar los conceptos de la geometría afín para diseñar un algoritmo, en lugar de complicarnos en 
desarrollos trigonométricos. Acudiremos al concepto de espacio vectorial 3D y haremos uso casi exclusivo de 
los conceptos de vector libre, producto escalar de dos vectores y producto vectorial de dos vectores. 


En lugar de una fórmula extensa y complicada intentaremos obtener un método por pasos, un algoritmo, que 
nos resulte cómodo para implementar en una aplicación móvil o una página web. 


Comenzamos enfocando nuestro objetivo: queremos 
hallar el rumbo inicial, azimut o ángulo que forma una 
línea imaginaria en el suelo que pisamos dirigida hacia 
el norte y otra línea imaginaria en el suelo que apunta 
en la dirección de la trayectoria que seguimos. 


La figura 5 ilustra lo que buscamos. Las dos flechas 
negras representan dos vectores sobre el plano del 
terreno. 


Uno de esos vectores es la proyección sobre el terreno 
de un vector que apunta hacia el polo norte. El otro 
vector es la proyección sobre el terreno de nuestra 
trayectoria de salida. 


Mirando de nuevo la figura 4, identificamos tres planos 


que nos resultan especialmente interesantes. 


dirección norte 


proyectada 

en el suelo 
dirección 

en que 
avanzamos 
punto que 
pisamos 
azimut 


El primero de ellos es el plano del terreno cuando pisamos el punto A, que resulta normal (perpendicular) al 
radio terrestre que va desde el centro O de la esfera hasta el punto A de partida. 


El segundo plano interesante es el plano vertical que contiene al punto A inicial y al eje z. 


El tercer plano considerado es el que incide a la vez con el punto A inicial, el punto B final y el centro O de la 


esfera. 


La intersección del primero y segundo planos origina una recta que contiene al vector que en el suelo apunta 
hacia la dirección norte. Y la intersección entre el primero y tercero planos engendra una recta que contiene al 


vector que sobre el suelo sigue la trayectoria inicial. 


Consiguiendo estos dos vectores calcularemos el ángulo que extienden y éste será nuestro azimut buscado que 
determinará nuestro rumbo inicial en la trayectoria del punto A inicial hasta el punto B final. 


En lugar de plantear sistemas de ecuaciones que expresen estas intersecciones entre planos, optaremos por 
construir vectores que apunten en las direcciones que nos interesan. Para ello vamos a dotarnos de unas 


funciones básicas que realicen las siguentes operaciones: 


página 9 de 24 


Fórmulas simples para navegación ortodrómica O 2022 Carlos Grasa Lambea 


Calcular el producto escalar de dos vectores, calcular el producto vectorial de dos vectores, calcular el vector 
unitario paralelo de otro y calcular el ángulo entre dos vectores. 


Si tenemos dos vectores genéricos, podemos nombrar al primero como uu mientras que el segundo será V 
en nuestra notación. Cada vector tiene tres componentes, sus proyecciones sobre los ejes x , y , Zz que 


en caso del vector Uu serán Xy , Yy , Z( como análogamente el vector Y tendrá sus coordenadas 
cartesianas Xy , Yy , Zy Correspondientes. 


Usaremos el símbolo + para el producto escalar de dos vectores, el mismo que la multiplicación, que será 
reconocido porque opera sobre vectores. Ejemplo: uv 


El símbolo A representará el producto vectorial de dos vectores. Así: UAV 
A la longitud de un vector le llamamos módulo del vector. El vector ú tiene su módulo Jú| asociado. 


> 


Por repasar conceptos de geometría, diremos que el resultado del producto escalar de dos vectores u-v esun 
número escalar equivalente al producto de la longitud (módulo) de uno de ellos por la proyección del otro sobre 
el primero. Así, podemos calcular con 

fórmula 17: — ú-v=|úl-[v|-cosw 

donde «w  eselángulo que forma ú con Y enel espacio. 

Si usamos las componentes cartesianas de esos vectores, hay otra forma de calcular su producto escalar con 
fórmula 18: — UÚV=Xy"Xy*+Yy'Yy+Zy"Zy 


Miremos qué ocurre al calcular el producto escalar de un vector consigo mismo: 


Ú=Xy"Xyt Yu YytZytZy=Xy + yv + zy=1ú|Jú)-cos0=|ú"=ú* 


sy 


Esto nos permite calcular el módulo de un vector mediante la 
fórmula 19: lúl=V xo + yo + Zo 

pero también mediante 

fórmula 20: la[=Vuú-u 


Esto va a permitirnos calcular un vector de longitud unidad que sea paralelo a otro: 


fórmula 21: 1,= 


Donde 1, significa un vector de longitud unidad paralelo a ú y esta fórmula nos permite reducir un 
vector para que tenga longitud 1, un vector unitario, tomando sus coordenadas cartesianas divididas entre la raíz 
cuadrada del producto escalar consigo mismo. 


Si ahora combinamos la 


fórmula 17: — U-v=|ú||v|-cosow 
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junto con la 

fórmula 20: lúl=Vú-ú 

sustituyendo y simplificando obtenemos: 
ú-v=vVú-ú-/v-V-cosw 

Ahora despejamos el coseno: 


cosu= 


ey] y 


Y ú 
vv V(u-0)(v-v) 


Vis 
Y llegamos a la ecuación que nos permitirá obtener el azimut, nuestro rumbo inicial, mediante la 


fórmula 21: co=arccos ( 


Además de profundizar en el producto escalar también vamos a estudiar otra operación entre vectores que 
usaremos en nuestros cálculos: el producto vectorial. El símbolo con que lo denotaremos será A que también 
se emplea en lógica booleana. 


El resultado de aplicar el producto vectorial a dos vectores arroja un tercer vector que es perpendicular (o 
normal) a los dos primeros vectores a la vez. El módulo del vector resultante es fácil de conocer mediante la 


fórmula 22: 4 v|=|úl-|v| senw 


Es similar a la fórmula del producto escalar, pero en este caso debemos calcular explícitamente el módulo y 
usamos la función seno. 


En el caso de usar las coordenadas cartesianas de los vectores también tenemos esta otra 


> > > 


fórmula 23: ÚAV=[ y y:Zy—Zy Y y) 1, +(Z "Xy—Xy"Z y) 1,+(Uy Vy—Uy Vx): 1, 


> 


donde el resultado es otro vector y denotamos 1, , 1, , 1, como vectores unitarios con la dirección de 
su respectivo eje. También podemos expresar el vector separando sus coordenadas mediante comas y 


encerrando la terna entre paréntesis, como en la 


fórmula 24: ÚAV=( Y y Zy—Zu'Yv ) Zu'Xy—Xy"Zy ) Ux"Vy—Uy'Vx) 


Hay que recordar que la dirección del vector resultante es 
siempre perpendicular a los dos vectores de partida. 


El sentido de la dirección viene determinada por la regla del 
tornillo dextrógiro (a derechas) según la figura 6. 
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Si solo nos interesase el módulo del vector resultante, podríamos aplicar, por ejemplo: 


fórmula 25: ¡BAD =V (Y yzy—2y Y y) +(Zy'Xy— Xy 2y YH Xy" Y y — Y y Xy) 


Trayendo de nuevo la 
fórmula 22: lA v|=|ú|-v]| sen w 
Podemos despejar el ángulo así: 


[A v]| 
july 


) 


fórmula 26:  «w=arcsen( 


Una vez provistos con estas herramientas, vamos a analizar 


cómo conseguir los dos vectores tangentes al terreno de la A 
figura 5 que determinarán el azimut que buscamos. tias 
en el suelo 


Al vector que apunta hacia el norte lo nombraremos Ññ ció 
(recordando norte), mientras que para el vector que marca la en que 

. ., A: r > avanzamos 
dirección del rumbo inicial usaremos el símbolo  r, 


(recordando rumbo inicial). 


Construiremos nuestros vectores sin emplear cálculos que usen 

divisiones con posibles problemas de infinitos. En vez de ello, 
aprovecharemos la propiedad del producto escalar para 

conseguir vectores unitarios, y la propiedad del producto 

vectorial para conseguir vectores normales (perpendiculares) a algunos planos. 


Vamos a identificar en la figura 7 algunos elementos 


que usaremos en adelante. polo 
norte 


Tenemos la esfera con su centro (O conectado 
mediante el vector ñ con el polo norte. 


Vemos representados el punto A - inicial y el punto 
B final, ambos conectados con el origen 


mediante los vectores d , b respectivos. 


Vamos a resaltar parte de un plano que contiene a los 


vectores dá , b como también el sector circular 
definido por los tres puntos O , A ,polo norte. 


punto que 
pisamos 


azimut 
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La zona gris de la figura 8 muestra el plano al que nos estamos polo 


norte 


refiriendo. 

Este plano es interesante porque pasa por el eje z que incide en el 

polo norte. Además es normal al plano del terreno en el punto A 

inicial. ¡ 
Un objetivo nuestro es encontrar un vector que, partiendo del punto 

se dirija hacia el norte y sea paralelo al plano del terreno. A este 

vector le llamaremos vV, pensando en que será el 0 “N del azimut. 

En la figura 9 hemos añadido el plano del terreno conteniendo al 


vector Vy (en color naranja) saliendo del punto A dirigido 
hacia el norte local del terreno. 


También hemos representado un vector v, (en color verde) que es perpendicular al plano gris, por lo cual 
también es normal a los vectores d , mM  intuyendo que podemos construirlo mediante el producto vectorial 
de estos últimos. 


Hemos ampliado el detalle en la figura 9b para mayor claridad. 

Vemos, entonces que el vector lo podemos obtener mediante la 

fórmula 27: — v,=ñAd 

Las coordenadas de ñ son las del polo norte: (0,0,R ) siendo R el radio terrestre. Pero nosotros 
buscamos como resultado un ángulo que no dependerá del radio de la esfera. En algún momento R desaparecerá 
en alguna simplificación. Por tanto, vamos a suponer el radio R como 1 sabiendo que esto no influirá en el 


resultado y que hace más facil el desarrollo. 


Es decir: 
fórmula 28: n=(0,0,1) 


Y recordando lo ya visto en la sección 1 de este documento: 


fórmula 29: á= (cosó,:cosA, ,cosó,send, , send, ) 
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También observamos que el vector v, es normal a los vectores 


fórmula 30: — Vy=dAVvV, 


que en combinación con la 


fórmula 27: — v,=ñAd 
llegamos a la 
fórmula 31: N=4M(ñAA) 


Ahora nos vamos a ocupar de conseguir un vector paralelo al terreno 


que señale la dirección en que avanzamos. 


d , V, expresándolo así: 


dirección norte 
proyectada 
en el suelo 


dirección 
en que 
avanzamos 


punto que 
pisamos 


azimut 


Este vector nos dirige por la ruta más corta hacia el punto B- final. Le asignamos la notación V, 


ilustrándolo en la figura 10. 


De manera similar a cómo hemos actuado en el caso de v, usaremos el producto vectorial para conseguir un 


vector como V¿ e€nunos pocos pasos. 


Para ello nos fijamos en la figura 10. Queremos 
conseguir un vector V¿ mostrado en color rojo, 


coplanario con el terreno donde también está vVy 


mostrado en naranja. El rumbo inicial estará 
determinado por el ángulo que formen ambos 
vectores. 


En primer lugar construiremos el vector  v, 


perpendicular al plano (en gris) incidente sobre los 
puntos A , O , B como apreciamos en la 
figura 10. 


Este vector V, será normal tanto al vector d 


como al vector b  intuyendo que el producto 
vectorial de estos últimos producirá el primero: 


fórmula 32: v,=4nb 


Y en un segundo paso contruimos el vector V, que 
recurriendo de nuevo al producto vectorial de estos 


es normal tanto al vector V, como al vector d 
últimos: 
fórmula 33: — V¿=V,Aá 


Sustituyendo en la ecuación anterior llegamos a la 


fórmula 34: va=(ánb)á 


Consiguiendo nuestro objetivo de dos vectores que extienden nuestro azimut buscado. 
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Por último vamos a exponer los pasos de nuestro algoritmo para hallar el rumbo inicial de la trayectoria entre 
los puntos inicial y final: 


Tendremos que solicitar al usuario las coordenadas geográficas de los puntos A, B como también asignar las 
coordenadas del polo norte: 


ñ=(0,0,1 ) aplicando la fórmula 28 


á= (cosó,'cosA, ,cosó, send, , send, ) aplicando la fórmula 29 


Usaremos una variable para cada coordenada en nuestro algoritmo. También consideraremos que las funciones 
trigonométricas en los lenguajes de programación usan radianes y convertiremos grados sexagesimales a 
radianes usando la fórmula: 


_ 


di = grados: 
radianes= grados: y 


Para terminar el algoritmo tenemos que calcular el ángulo que forma el vector v,  conel vector vz 
pudiendo optar entre varias posibilidades. Por un lado podríamos usar la ya vista 


fórmula 21: c=arccos ( 
aplicándola a los vectores V, , Vy teniendo en cuenta que son vectores unitarios y, por tanto, ya no es 
necesario calcular sus módulos para el denominador. 


Por otro lado podemos calcular el seno y el coseno del azimut para hacerlo de una manera con menos 
ambigúedades, como mostramos a continuación: 
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Calculamos el coseno mediante la ya conocida 

fórmula 17: U-v=|ú||v|-cosow 

sabiendo que los vectores Vy , VÉ¿ sonunitarios y simplificando adaptamos la expresión: CosWw=Vy'Vy 
Como trabajamos con las cordenadas cartesianas, hacemos uso de la 

fórmula 18: — UÚ-V=xX¡¿"Xy+Yy Y y+Zy"Zy 

adaptándola a nuestros vectores, formando la 

fórmula 35: OSI=A AAVV FL 

De manera similar obtendremos el seno del azimut, esta vez usando el producto vectorial, para ello recordamos: 


fórmula 22: lún v|=|úl-|v|senw 


fórmula 25: ¡UA v|= Vw Zy=Zy Y y) HZ ycXy—=Xy 2 HU Vy Uy V y) 


Las combinamos y adaptamos a nuestros vectores, teniendo en cuenta de nuevo que son unitarios para formar la 


fórmula 36: seno=W[YwZa—Zw Ya HZ XX Za) + Y p— Y vv X a) 


En la figura 11 ilustramos el azimut en color verde junto con los 
vectores que lo engendran. 


La mayoría de lenguajes de programación disponen de una 
función atan2() que es muy útil para pasar de coordenadas 
cartesianas a polares. 


Esta función requiere dos argumentos: el seno y el coseno, dicho 
de otro modo las coordenadas cartesianas x , y cuando queremos 
obtener el equivalente a 


arctg Ln 


La notación, por ejemplo, en Javascript de esta 
función es: 


fórmula 37: azimut=atan 2( yx) 


Y con esto ya podemos plasmar nuestro algoritmo ordenado por pasos y comentado. 
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Algoritmo para calcular el rumbo inicial de una trayectoria ortodrómica entre dos puntos conociendo sus 
coordenadas geográficas: 


paso 1: introducir latitud del punto A en: 6, 
paso 2: introducir longitud del punto Aen:A, 
paso 3: introducir latitud del punto Ben:6 y 
paso 4: introducir longitud del punto B en: Az 
paso 5: o 

paso 6: A=Am 0 

paso 7: 0=9»380 

paso 8: Ag= a 

paso 9: Xy=0 

paso 10: yy=0 

paso 11: Zy=1 

paso 12: X¿=C080,'C0SA, 

paso 13: y ,=c0sÓ send, 

paso 14: z =senó, 

paso 15: Xg=C0SÓ ¿'COS Ay 

paso 16: y¿=C08Ó ¿'SenA y 

paso 17: z,=senóz 

paso 18: AI IV 

paso 19: Val X= XX Za 

paso 20: OOO ES 

paso 21: Xy Y a Ly1 Za Y y 

paso 22: Va Za A 
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paso 23: Zi XA Y y Y a Xy 
MARA 
paso 24: Mn= lt +2 
Xvy 
paso 25: X vw ==— 
VN Pr 
Y vn 
aso 26: JERER 
p Y vn My 
Zn 
paso 27: Ziy= 
VN Mo 
paso 28: Xv2= YA Zg ZA YB 
paso 29: Yy2224 Xp Xp Zp 
paso 30: Zvo=XA YB Ya Xp 
paso 31: XvB=Yvy2 ZA" Zly Ya 
paso 32: Yva=2Zy2XA4 7" %Xy2ZA 
paso 33: ZiuB=*Xv Ya Yy2Xa 
EA 
paso 34: Mys= x2,+ Va +tZm 
Xvyg 
paso 35: A 
VB EN 
Y va 
36: == 
paso Y va a 
Zvg 
paso 37: Lia 
mw 
paso 38: 0x=X yw Xyg* Y yy" Y yg +Z yn Z yg 
paso 39: =V (Y acZys—Zyw Y vo) +(Zy X ya X Za) Hr Y vs Y var Ava) 
paso 40: o=atan2(o, o) A 
paso 41: Escribir rumbo inicial R, desde: w 
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Sección 3: Comprobación de la fórmula de la distancia ortodrómica y del algoritmo del rumbo inicial. 


En esta última sección vamos a poner a prueba nuestros desarrollos de la fórmula para el cálculo de la distancia 
ortodrómica y del algoritmo para el cálculo del rumbo inicial de la trayectoria ortodrómica. 


Vamos a tomar las coordenadas geodésicas de dos ciudades, Zaragoza en España y Berlín en Alemania, para 
luego comprobar nuestros resultados con los de alguna célebre calculadora online pública para navegación: 


Zaragoza 
latitud: 41.65078071020651 “N 
longitud: -0.8888014436201552 “E 
Berlín: 
latitud: 52.520779305747965 “N 


longitud: 13.38960953926479 "E 


Radio terrestre: 6371000 metros 


Asistidos por una hoja de cálculo realizamos el cálculo de la distancia ortodrómica entre las ciudades escogidas 
y con el radio terrestre citado. Introduciendo los datos de coordenadas geográficas y las fórmulas de conversión 


a radianes y de cálculo de la distancia. 


El resultado aparece en la figura 12: 


E | latitud: | 41,6507807102065 
9 | longitud: | -0,888801443620155 BE 


AN] 
[Er 
Latitud: | 52.520779308748 | *N>>] 
—hongitud: | 13.3896095392648 
metros | 


6371000 
distancia ortodrómica: 1615757,61955421 | metros | 


0,726943259419246 
-0,0155125115876505 


0,916660524598604 
0,233692772016612 


Vamos a comparar este resultado de 1615757.61955421 con el de alguna página web pública y fiable, como 


http://rodamedia.com/navastro/online/javascripts/ortodromica_comp.htm 


En la figura 13 mostramos un pantallazo de una parte de la página web: 
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¡A] O A rodamedia.com/navastro/online/javascripts/ortodromica_comp.htm 200% YY Q Buscar XxXbBa QM DO a4s= 


Navegación ortodrómica 
Origen (Situación inicial): 


Latitud: — [41.650 grados [0 minutos — [Norte |[N/S] 
Longitud : [0.888€ grados |0 minutos — [Oeste EN 


Destino (Situación final): 


Buenas noches 


Miércoles, 27 /07/2022 


3:34:36 UTC Latitud:  [52.52C|grados [0 minutos — |Norte M3 
Longitud: |13.38S| grados lO | minutos Este (EW) 


Inicio l Calcular || Borrar | [Volver | 


Publicaciones 
Navegación 
Meteor: logía 


Resultados: 


Astronomía 
Rumbo ortodrómico inicial: 136.735679 | grados 
Distancia ortodrómica: 1871.8514423: millas náuticas 


Ahora convertimos la distancia ortodrómica de millas a metros para obtener el resultado. También tomamos 
nota del rumbo inicial ofrecido en la misma página web. 


distancia ortodrómica: 1615757.61955421 metros 
rumbo ortodrómico inicial: 36.7356798707116 *“N 


Observamos que el cálculo de la distancia coincide con el que nosotros hemos realizado perfectamente. Ello es 
porque están empleando la misma expresión que nosotros: 


fórmula 15:  d=R-arccos(senó,:sen$,+c0só,:cos 9ycos(14—Ag)) 


Acudimos a una segunda página web pública y de confianza, que 


A z z a . > . Ortodrómica - Great Circle Saili: 
realiza un cálculo más sofisticado, no en el cliente con Javascript sino di 


en su servidor Acepta escribir: grados, minutos y decimas de minutos (dejar en blanco segundos) o 
grados, minutos, segundos y decimas de segundos 
https://www.webmar.com/web/calc/cal_orto.php E 
Latitud : 
41.65 grados O minutos O segundos Norte |N/S 


Longitud : 
0.881 grados O minutos lO segundos Oeste ||E/W| 


En la figura 14 vemos parte de la página web donde hemos introducido 


los mismos datos y como resultado obtenemos: ds 


Latitud : 


distancia ortodrómica: 871.8514423157484 millas náuticas sae groso Immaosio segundos None NS 

rumbo inicial: 36.735679870711564 “N A A 
E Calcular Restablecer 

Convirtiendo las millas náuticas a metros con solo multiplicar por los 

1852 metros que tiene cada milla náutica, llegamos a una iaias 

distancia ortodrómica: 1614668.8711687660368 metros 5075507] vrados 


Distancia Ortodromica: 


871.8514423 millas nauticas 
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y con nuestro resultado: 1615757.61955421 metros 
estimamos nuestro error relativo: 0.0674% 


Este error relativo nos parece muy aceptable para nuestros propósitos. 


Las dos páginas web arrojan un resultado prácticamente idéntico para el rumbo inicial; nos quedaremos con la 
referencia de 36.735679870711564 “N para comparar con el rumbo inicial que obtengamos nosotros. 


Vamos a codificar nuestro algoritmo en lenguaje Javascript dentro de un documento HTML mínimo 
añadiéndole el cálculo de la distancia ortodrómica. 


Programa en Javascript para calcular la derrota y el rumbo inicial de una trayectoria ortodrómica entre dos 
puntos conociendo sus coordenadas geográficas: 


<!DOCTYPE html> 
<html lang=es > 


<head><title> Cálculo del rumbo inicial ortodrómico</title> 


<meta charset="utf-8" /><meta name=viewport content="width=device-width, initial- 
scale=1.0" /> 


<meta name=author content="0 2022 Carlos Grasa Lambea" /> 


<style> 
.separal ( 
display: inline-block; 
width: l1l5rem; 
text-align: right; 


) 
</style> 


</head> 

<body> 
<h1>Cálculo de distancia y rumbo inicial ortodrómicos</h1> 
<hr /><br/> 


<span class=separal >latitud A: £nbsp;</span><input id=latituda type=text size=25 
value=41.6568 /><br/> 
<span class=separal >longitud A: £nbsp</span><input id=longituda type=text size=25 
value=-0.8785 /> 

<br/><br/> 

<span class=separal >latitud B:énbsp</span><input id=latitudb type=text size=25 
value=52.5163 /><br/> 
<span class=separal >longitud B:s8snbsp</span><input id=longitudb type=text size=25 
value=13.3777 /> 

<br/><br/> 


<span class=separal ></span> 

<input id=botoncalcular type=button value="calcular" onclick="calcular ();" /> 

<br /><br/> 

<span class=separal >distancia ortodrómica:s£nbsp;</span><input 
id=distanciaortodromica type=text size=25 value="" /><br/> 


<span class=separal >rumbo inicial ortodrómico:gnbsp;</span><input id=azimut 
type=text size=25 value="" /> 


página 21 de 24 


Fórmulas simples para navegación ortodrómica O 2022 Carlos Grasa Lambea 


<br/><br/><hr/><br/> 


<section id=1lo0gs ></section> 
<script> 
// Algoritmo para cálculo de Rumbo Inicial Ortodrómico mediante productos escalar y 
vectorial 
// 0 2022 Carlos Grasa Lambea 
const calcular=function() ( 


let r=6371000; // radio volumétrico medio terrestre en metros 


// introducción de datos comprobando validez 


let latA=document.querySelector ("tlatituda") .value; // latitud del punto inicial 
let lonA=document.querySelector ("+longituda") .value; // longitud del punto inicial 
let latB=document.querySelector ("tlatitudb").value; // latitud del punto final 
let lonB=document.querySelector ("+longitudb").value; // longitud del punto final 
// conversión de grados sexagesimales a radianes 


latA*=Math.P1/180; 

lonA*=Math.P1/180; 

latB*=Math.P1/180; 

lonB*=Math.P1/180; 

let n=[0, O, 1]; // coordenadas del polo nort n esfera de radio 1 
let a=cartesianas (latA, lonA); // vector del punto A en coordenadas cartesianas 
let b=cartesianas (latB, lonB); // vector del punto B en coordenadas cartesianas 
let vli=productoVectorial (n, a); // cálculo del vector vl como producto vectorial de 
na 

let vn=productoVectorial (a, vl); // vector vn como producto vectorial de a”"vl 
vn=unitario (vn); // hacer unitario el vector vn 

let v2=productoVectorial (a, b); // vector v2 como producto vectorial de a”b 

let vb=productoVectorial (v2, a); // vector vb como producto vectorial de v2%a 
vb=unitario (vb); // hacer unitario el vector vb 


// calcular coseno y seno del azimut 


let wx=productoEscalar (vn, vb); // coseno como producto escalar de <vn,vb> 
let wy=modulo (productoVectorial (vn, vb)); // seno como módulo del producto vectorial 
de vn*vb 


let rio=Math.atan2 (wy, wx)*180/Math.Pl; // azimut del rumbo inicial ortodrómico en 
grados 
if (rio<0) // si es negativo, añadir 360 

rio+=360; 
let dort=r*Math.acos (Math.sin(latA) *Math.sin(latB) 
+Math.cos (latA) *Math.cos (1latB) *Math.cos (1lonA-1onB)); // distancia ortodrómica en 
metros 


// imprimir distancia ortodrómica y rumbo inicial ortodrómico 
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document .queryS 


metros"; 


lector ("Fdistanciaortodromica") .value=Math.round (dort*10000)/10000+" 


lector ("fazimut") .value=Math.round (rio0o*10000)/10000+" *N"; 


document .queryS 


return; 


y; 


const productoVectorial=function(a, b) ( 
return [a[1]*b[21-a[2]*b[1], a[l2]*b[01-a[0]*b[2], a[0]*b[11-a[1]*b[0]]; 


l; 


const cartesianas=function(lat, lon) ( 
return [Math.cos(lat)*Math.cos(lon), Math.cos (lat) *Math.sin(lon), Math.sin(lat)]; 


H; 


const unitario=function(v) ( 


let norma=modulo (v); 
if (Math.abs (norma)<le-9) 


return v; 


else 


return escalarVector (1/norma, v); 


y; 


const productoEscalar=function(a, b) ( 
return a[0]*b[0]+a[1]*b[1]+a[2]*b[2]; 


H; 


const modulo=function (v) ( 


return Math.sqrt (productoEscalar (v, v)); 


; 


const escalarVector=function(k, v) ( 
return [k*v[0], k*v[1], k*v[2]]1; 


y; 


const sumaVectorial=function(a, b) ( 
return [a[0]+b[0]1, a[1]1+b[1], a[2]+b[2]1; 


y; 


const restaVectorial=functionla, b) ( 


return sumaVectorial (a, escalarVector (-1, b)); 


y; 


</script> 


</body> 
</html> 
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Fórmulas simples para navegación ortodrómica O 2022 Carlos Grasa Lambea 


El programa toma por omisión las coordenadas geográficas de la ciudad de Zaragoza en España como punto 
inicial y las de Berlín en Alemania como punto final. 


La ejecución del programa nos devuelve el siguiente resultado 


Distancia ortodrómica: 1615757.619554206 metros 
Rumbo inicial ortodrómico: 36.735679870711564 “N 


Lo cual coincide perfectamente con las referencias tomadas de las páginas web consultadas. Vemos así que 
nuestro método es aceptable para las aplicaciones que tenemos en expectativa. 


Carlos Grasa Lambea 
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